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Abstract 

Among the multiple steps constituting the kinesin's mechanochemical cycle, one of the 
most interesting events is observed when kinesins move an 8-nm step from one microtubule 
(MT)-binding site to another. The stepping motion that occurs within a relatively short 
time scale (~ 100 ^s) is, however, beyond the resolution of current experiments, therefore 
a basic understanding to the real-time dynamics within the 8-nm step is still lacking. For 
instance, the rate of power stroke (or conformational change), that leads to the undocked- 
to-docked transition of neck-linker, is not known, and the existence of a substep during the 
8-nm step still remains a controversial issue in the kinesin community. By using explicit 
structures of the kinesin dimer and the MT consisting of 13 protofilaments (PFs), we 
study the stepping dynamics with varying rates of power stroke (k p ). We estimate that 
kp 1 < 20 fis to avoid a substep in an averaged time trace. For a slow power stroke 
with kp 1 > 20 fj,s, the averaged time trace shows a substep that implies the existence 
of a transient intermediate, which is reminiscent of a recent single molecule experiment 
at high resolution. We identify the intermediate as a conformation in which the tethered 
head is trapped in the sideway binding site of the neighboring PF. We also find a partial 
unfolding (cracking) of the binding motifs occurring at the transition state ensemble along 
the pathways prior to binding between the kinesin and MT. 

Introduction 

First recognized via their close relationship between ATPase activity and organelle transport 
along MTs 0], kinesins have received broad attention as a prototype of molecular motors 
for the last two decades. Recent single molecule (SM) experiments have shown that, with each 
stepping motion being strongly coupled to the ATP, the kinesin moves toward the (+)-ends 
of MTs by taking discrete 8-nm steps in a hand-over-hand fashion [7j, a a. 

Although the ultimate understanding to the kinesin's motility is still far from the completion, 
the SM experiments together with the series of kinetic ensemble measurements 



id . I 111 . 1 121 . 1 1 31 ] and theoretical studies , [l5|, H, [13, H, [l9| begin providing glimpses to the 
physical principle of how kinesin walks. 

Along the kinesin's mechanochemical cycle (SI Fig. 7), one of the main observations of 
the SM experiments is the stepping dynamics that enables the kinesin to move forward. The 
actual time spent for the stepping motion itself (< 100 /xs), compared to the ATP binding and 
hydrolysis (> 10 ms), is however too short to detect the details of the dynamics with the spatial 
and temporal resolution of current instruments. Thus, it is still difficult to answer many basic 
questions [20J related to the stepping dynamics. Some of those questions are : (a) During the 
8-nm step, how does the swiveling motion of the tethered head occur? Does any detectable 
substep exist that reflects a transient intermediate? 21, 22, 0] (b) What fraction of the time 



and length scales is contributed from the power stroke and the diffusional search? (c)_ Does the 

To shed 



24 



kinesin walk parallel to the single PF or walk astride using two parallel PFs? [23, 
light on these questions, we propose to take advantage of the native topology of kinesins and 
MTs. For instance, our previous study has clarified the regulation mechanism between the two 
heads [25|, |57J] by using the native topology-based two-head bound model of kinesin on the MT 



19(| . Following the same line of thought, we show that even the dynamical pathways of kinesins 



reflect "the topological constraints emanating from the molecular architecture." In the present 

to study the stepping dynamics 
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work, we have adapted our previous two-head bound model 
of kinesins on the 13-PF MT. 

According to the "neck-linker docking model 26|, |27|, |28|," the stepping motion of kinesin is 
initiated by the undocked-to-docked transition of the neck-linker at the MT-bound head (power 
stroke), and the rest of the binding process of the tethered head to the next MT-binding site is 
accomplished via a diffusional search. The whole stepping dynamics is typically interpreted as a 
combination between the directional motion of power stroke and the nondirectional diffusional 
search. At the molecular level, the length of the neck-linker shrinks gradually, with a progress 
of power stroke, making transition from a disordered to an ordered state; the resulting configu- 
rations bias the diffusional motion of the tethered head in one direction. Thus, we can rather 
interpret the power stroke as a time dependent boundary condition that biases the diffusional 
motion, considering that the tethered head diffuses on the time dependent energy landscape to 
search for the next MT-binding site. 

To study the dynamics, we first obtain through simulations potentials of mean force (PMF) 
felt by the kinesin's tethered head on the MT at two extreme cases; one is the PMF with 
fully disordered neck-linker (F\ =0 (x, y, z)) and the other with ordered neck-linker (F\ = i(x, y,z)), 
where A is a parameter that specifies the degree of the neck-linker being zippered. The kinesin 
structures with a disordered (A = 0) and an ordered neck-linker (A = 1) at the MT-bound 
head are generated by switching on and off the native contacts inside the green circles in Fig. 
IB (see Results and Methods for more details). We subsequently perform Brownian dynamics 
simulations of a quasi-particle representing the tethered head on a PMF that is mixed between 
the two extreme cases. The disordered-to-ordered transitions are mimicked by mixing the two 
PMFs, F\ =Q (x, y,z) and F\ = i(x, y, z), in a time-dependent manner. We show that because of 
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the multiple binding sites on the MT surface, the tethered head has a chance of misbinding to 
other MT-binding sites. To avoid such an intermediate, we argue that the rate of power stroke 
(k p ) should be faster than the sampling rate on the MT surface Throughout the paper, 

we designate the rate of power stroke and the rate of space exploration by diffusion as k p and 
fee, respectively. 

The present work provides scenarios on how the dynamic pathways of the tethered head 
depends on the rate of the power stroke, based on the landscape of stepping and the molecular 
details that the current experiment cannot easily access. Importantly, the present work will give 
further insights to resolve the recent experimental debate |22l . 
substep within the kinesin's 8-nm step. 
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20 on the existence of a 



Results 

Topological constraints in stepping dynamics. Upon ATP binding to the empty MT- 
bound head (positioned at the binding site o in Fig. 1A), the allosteric communication between 
the strained nucleotide binding pocket and the disordered neck-linker leads to the docking of 
the neck-linker to the (37 and L10 neck-linker binding motifs. As a result, the tethered head 
swings forward from one MT-binding site (a in Fig. 1A) to the other (e in Fig. 1A) while 
the MT-bound head remains at the same position (see SI Fig. 8 for the snapshots of a kinesin 
during the stepping). For the tethered head to bind to the MT-binding site, the ruggedly shaped 
kinesin head having specific MT-binding motifs should explore the MT surface and fit into the 
MT-binding site in a right orientation under topological constraints. The two neck-linkers (one 
from the MT-bound head, and the other from the tethered head) and the steric hindrance be- 
tween the two heads restrain the search space available for the tethered head. The topology of 
MT also makes unique the association process between the tethered head and the MT-binding 



site. In the MT made of 13 PFs |29J the lateral binding interfaces between the adjacent PFs are 
shifted by ~ 0.95 nm, which makes the MT helical (Fig. 1A). The distances between the diag- 
onally located tubulins become asymmetric (|o<i|(= 9.4 nm) < \of\(= 10.9 nm)). The distance 
between the interfaced tubulins at the side (\6c\ = \dg\ = 6.3 nm) is shorter than the distance 
between the tubulins along the same PF (|de| = 8.0 nm). In light of the contour length of a 
neck-linker (15 aa x 0.38 mxv/aa ~ 5.5 nm), two fully extended neck linkers (~ 11 nm) allow the 
tethered head to cover a wide range of the MT surface. Since every MT-binding site on the MT 
equally interacts with the tethered head, one cannot totally rule out the possibility of binding 
to the neighboring PFs while the experiments suggest that the kinesin move straight along the 
MT 
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PMF between kinesin tethered head and MT. During the disordered-to-ordered tran- 
sition, the length of disordered neck-linker decreases gradually. Depending on the length of 
the disordered neck-linker, there are substantial variations in the search space available for the 
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tethered head (Fig. 2A and Fig. 3A). We sample the conformational space using the centroid 
of tethered head from the multiple trajectories with varying temperatures (see Methods for the 
energy function), and calculate the 2-D PMF between the tethered head and the MT, projected 
on the xy, xz-, yz-planes, at two extreme cases (see SI text for the computational details of PMF 
construction). One results from search processes under which the neck-linker of MT-bound head 
is practically in the zippered state (A = 1). We expect that this process is realized when the 
neck-linker zippering rate is faster than the rate of space exploration (k p ^> Whereas, the 
other results from the case when the neck-linker is in disordered state throughout the search 
process (k p <^ fee, A = 0). Although the exact values of k p and &g are not known, we presume 
that the stepping dynamics should occur on a PMF linking between these two regimes. 

From the structures with A = 1, we find two major basins of attraction in the 3-D space (Fig. 
2). One is the MT-binding site e and the other is a broad basin (S) formed at the forward-right 
corner relative to the MT-bound head (Fig. 2C, Fi(x,z)). Between the two basins exists the 
free energy barrier of ~ 1 — 2 ksT in Fi(x, z) or Fi(x, y), and ~ 4 — 5 ksT in Fi(y, z), depending 
on the projection. The basin S stems from the large conformational degrees of freedom that the 
tethered head should explore before reaching the binding site e (Fig. 2C, circles in blue). The 
energetic bias to the binding site d (Fig. 2C, arrows in magenta) is also present, but is negligible 
(< 1 ksT). With an increasing temperature the stability of basin S relative to the binding 
site e increases, which indicates that the nature of the broad basin is entropic (Fig. 2D). Note 
that the 1-D free energy profiles (Fig. 2B) projected from higher dimensional representation of 
free energy surface underestimates the barrier height by merging the various possible dynamic 
pathways. 

From the structures with A = 0, which is obtained by making repulsive the attractive neck- 
linker zipper contacts (the native contacts inside the green circles of Fig. IB, see also Methods), 
we find that the most frequently visited binding site is the site c that belongs to the adjacent 
PF (Fig. 3). Explicit analysis on the bound complex shows that the binding of the tethered 
head to the site c is not as complete as the one to the site e. Interestingly, the structure at c has 
only ~ 1/3 of the interfacial native contacts (see Fig. 4). We find that a strained neck-linker 
induces a significant distortion in the a6 helix, preventing a full binding. 

Modeling the stepping dynamics. The most straightforward strategy to monitor the 
real-time stepping dynamics of our model is to integrate the equation of motion for each coarse- 
grained unit of kinesin molecule in an overdamped environment. The inclusion of hydrodynamics 
to the simulations, which is essential to naturally retrieve a correct behavior of the translational 
diffusion for the tethered head (see SI text with SI Fig. 9), is however computationally too 
expensive for a model with more than 700 coarse-grained units. Thus, we take an alternative 
method to study the real-time dynamics by using the two PMFs obtained above. 

Conceptually, the stepping dynamics coupled to the power stroke is considered as a diffusional 
process with a moving reflecting boundary condition that restrains the search space, which is 
illustrated in Fig. 5A. The blue shade and the red dots depict the energetic cost for extending 
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the neck-linker and the multiple MT-binding sites, respectively. The progress of the neck-linker 
zippering is depicted by dotted lines with parameter A. The reflecting boundary moves from 
A = 0toA = lina finite rate. Depending on the rate of the variation from A = to A = 1 
as well as the motility of the kinesin head, the dynamic pathway of binding process changes 
because of the influence of the sideway binding site c. The rate of variation from A = to A = 1 
is defined as k p = r" 1 = dX/dt, i.e., A = k p t with < t < r p . 

We propose that the PMF time-dependently linking between the two extreme cases be mod- 
eled using the following ansatz adapting the dual-Go potential 
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( _ log \(l _ ±\ e -m(x,y,z) + ± e -l3F 1 (x,y,z)] (o<t< Tp ) n . 
{PF^x^z) (t>T p ) 

where the switching rate between F (x,y,z) and F\{x,y,z) is controlled by t p (— k^ 1 ). We 
perform Brownian dynamics simulations of a quasi-particle, representing the centroid of teth- 
ered head, on the "time-dependent" PMF, F(x,y, z,t), from the initial position of the head 
at (xi,yi,Zi) = (—6.0,4.0,-2.0) nm. We update the position at time t using f(t + At) = 
f(t) - D e J f VF{x, y, z, t)At/k B T + R(t) where R(t) is a vector of Gaussian random number sat- 
isfying (R) = and (R ■ R) = 6D e J^At. We choose D% = 2 fim 2 /s for the effective diffusion 
coefficient of the tethered head to reproduce a similar time scale as the one in the recent SM 



experiment by Yanagida and coworkers [21], where they monitored the stepping dynamics of 
kinesin in the time resolution of ~ 20 fis and suggested that the kinesin take substeps at the 
displacement of Ax ~ 4 nm. 

Simulating the dynamic trajectories by varying the t p , we count the number of trajectories 
trapped in the binding site c (Fig. 5B). The fraction of trajectories directly reaching the site 
e decreases almost exponentially as t p increases. We find that when r p < 20 /is, 90 % of the 
stepping dynamics occur without being trapped to the intermediate. The signature of such an 
intermediate, manifested as a "substep" in the averaged trace, is captured only when r p > 20 
fis (Fig. 5C). Note the similar itybet ween the patterns of averaged time traces by Yanagida and 



coworkers (see Figure 4a in Ref. [211]) and the present results plotted in Fig. 5C. Interestingly, the 
averaged trace for t p = 20 fis fits to the double exponential function (Fig. 5C) where the time 
scale of the fast phase corresponds to that of the power stroke, and the slow phase represents 
the rate process to overcome the barrier from broad basin to the target basin shown in F\(x, y, z). 

Partial unfolding (cracking) of structure facilitates the binding process. In the 

above sections, we studied the association between the tethered kinesin head and the MT-binding 
site, using the centroid of the tethered head. When we probe the kinesin conformation and the 



binding interface using native contacts, we observe partial unfolding 32| of the MT-binding 
motifs in the kinesin head at the transition state ensemble prior to the complete binding. The 
flexibility of the structure eases the binding process by reducing the entropic barrier to over- 
come, which is well known phenomena in protein-protein or protein-DNA association process 



33j, [34j, [35|. To quantify the degree of cracking in the kinesin structure during binding, we 
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use the fraction of native contacts for the MT-binding motifs of kinesin (Q p ) and the fraction 
of interfacial native contacts between the kinesin and the MT (Qint) (Fig- 6). In an exem- 
plary trajectory (Fig. 6A) starting from Q p pa 0.75, the MT-binding motifs are disrupted to 
Q p ^ 0.65. When the kinesin is completely bound to the MT, Q p pa 0.82(±0.08). By collecting 
the configurations and applying the histogram re-weighting technique (see SI text), we obtain 
the 2-D free energy surface F(Q p , Qi n t). Figs. 6C, 6D show that the Q p value of the kinesin-MT 
complex is greater than the value of the transition state ensemble. At higher temperature, the 
trend of cracking prior to the binding is more pronounced, showing a downward curvature in 
the pathways connecting the separated molecules and the complex. Interestingly, the 1-D free 
energy profile F(Qi nt ) (Fig. 6B) shows that the free energy barrier for the binding is ~ 6 &bT, 
which is higher than the ones measured as function of spatial reaction coordinates in the Fig. 
2. 



Discussions 



To study the dynamics of kinesin's stepping motion, we first made full use of the topological 
information available from the structures of kinesin and MT to build the PMF, and second 



considered the whole stepping dynamics as a "rectified diflusional motion" 36j, |37|] by envisioning 
the power stroke as a moving reflecting boundary condition for the diflusional motion of the 
tethered head. Although the level of our description is solely hinged on the native topology, 
lacking in the chemical details, several key issues on the stepping dynamics can be discussed in 
a semiquantitative fashion. 

(i) The PMF between the tethered head and the MT showed that the leftward diagonal 
stepping (of) is forbidden because of structural constraints (see Figs. 2 and 3). Also, the 
alternation between the sideway stepping (6c stepping) and the parallel stepping would generate 
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a helical path in the long run, which contradicts to the previous experimental findings 
Thus, the likelihood for using two parallel PFs is ruled out. An interesting finding in this study 
is that the intermediate structure trapped into the sideway binding site c has only ~ 30 % 
stability compared with the correctly bound complex at e. We expect that under the action 
of power stroke, the intermediate state found for k p <C k E becomes further destabilized and 
readily loses its binding with the site c. A direct comparison between F (x, y, z) and F\(x, y, z) 
at (x,y,z) pa (0,4,5) nm indicates that the binding site c is destabilized by > 5 ksT upon 
A = -> A = 1. 

(ii) With the speed limit of protein folding rate ~ 0(1) (fis)~ l 38|, |39|, |40] and the number of 
amino acids consisting of neck-linker N ~ (12 — 15), one can roughly estimate the rate of power 
stroke using the scaling relation for the folding rate of proteins with N (kp ~ k° F exp (— l.liV 1 / 2 ) 



38jj kUjj or kp ~ fc£exp(-0.36iV 2 / 3 ) with (k\ 



[IS 



42). The 



with (kp)" 1 ~ 0.4 lis 

agreement of kp X (~ 20 — 70 fis) with /c~ 1 (~ 20 /is) shows that the power stroke is closely 
connected with the conformational change of neck-linker whose activation barrier is estimated 
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as AG x /k B T{^ l.liV 1 / 2 or 0.36iV 2 / 3 ) -2-4. 

(iii) Given the diffusion constant (Dj[ ) and the approximate shape of PMF, one can estimate 
k B , the exploration rate over the PMF, by approximating the kinesin's motion as a Brownian 
motion in a harmonic potential (F(x) ~ 1/2 x kx 2 , i.e., Ornstein-Uhlenbeck process), whose 
conditional probability is solved as follows 43 . 



W(x,t\x ) 



2nk B T 
k ' 



-2t/r E - 



1 -1/2 



exp 



(x — xoe 



-t/r E \2 



2kgT I -i 



-2*/te ~ 



(2) 



where 
Since fc p 



: k B T/D e J f k. Whent > r £ , W(z,*|zo) -> P eg (ar), (x) e(? -> 0, and ((5x) 2 ) e(? -> k B T/k. 
0.02 k B T/nm 2 from the fit of 1-D -Fq(x) to a quadratic potential (Fig. 3C), and 



2 [im 2 /s, one gets = 25 /is. The lower bound for the rate of power stroke (k p ) is 
similar to the upper bound for the exploration rate of harmonic potential = t^ 1 ). If the 
PMF switching is too fast, the molecules sample only the subregion of the landscape, reflecting 
the signature of the far-from-equilibrium dynamics. 

(iv) The experimentally measured time traces for the rising phase, averaged over the different 

showed the signatures of intermediates. 
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stepping time scales, by Yanagida and coworkers 
The comparison between the data by Yanagida and coworkers' and our time traces generated 
over varying k p suggests that the k p distribute broadly, given by a distribution of rate constant 
g(k p ). Because of the molecular origin of the power stroke, associated with a complex energy 
landscape representing the molecular architecture, it is natural to speculate that the stepping 



441 ] . The multiple basins (the binding 



dynamics occurs via the kinetic partitioning mechanism 
sites c, d, and entropic basin S) in the PMF suggests that kinesins reach the target binding site e 
along multiple parallel pathways. Not dividing the time traces into the different classes like the 



way Yanagida and coworkers adopted 2l| , one should be able to fit the full ensemble average of 



the time traces to the finite number of multi-exponential function as (x)(t) = J °° dk p g{k p )e~ kpt 
451 ] . Although a visual signature of intermediate (a substep in the averaged time trace) is 
masked if the contribution from the fast kinetics (r p < 20 fis) is dominant, a careful statistical 
analysis of (x)(t) averaged over the entire time traces would give a glimpse to the g(k p ) that 
encodes the underlying energy landscape associated with the stepping dynamics. 



Methods 

To generate the configurations of kinesin on the MT, we performed hybrid Monte-Carlo/Molecular 
Dynamics simulations by adapting an ADP-complexed crystal structure of rat kinesin dimer 
(PDB code : 3kin) on the 13-PF MT structure (see SI text for further details). The 13-PF 
structure is built with multiple tubulin dimers fitted to the Downing and coworkers' 8-A reso- 



lution electron density map 29 
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Energy function. The energy function for the kinesin/MT system is modeled as Vtot = 
V K + V k .ut where V K {= Vg + Vg + V^ T ) is the energy function for the kinesin dimer and Vk-mt 
is the interaction between the kinesin and the MT. B and T denote the MT-bound head and 
the tethered head, respectively. V^ T describes the interaction between the B and T head. 



To design each kinesin monomer we use the self-organized polymer model 46|, [31 



1=1 ^ 

7V e -3 
i=l j=i+3 



(r M+1 - r° l+1 ) 2 
Rl 



i.i 



i=l 



Si 



a 



r i,i+2 



N e -3 N{ 
i=l j=i+3 



13 1 



(3) 



where £ =B or T. The first term models the chain connectivity with k = 100 kcal/ '(mol-A 2 ), 
R Q = 2 A. The second term decribes the attraction between the native contact pairs i and 
j. The native pairs are defined using Ay = 1 for the residue pairs within R c = 8 A at the 
native structure, and r°- is the corresponding distance. Eh(i,j) controls the strengths of native 
pairs. We choose Eh{i,j) = 1.8 kcal/mol for the residue % and j(= % + 3) within an a-helix, 
and Eh(i,j) = 1.2 kcal/mol for all other cases, e.g., /^-sheets, loops, and residue pairs between 
the secondary structural elements. The third term prevents the volume overlap between the 
residues i and i + 2. The last term is for the repulsive potential between the non-native pairs, 
modeled using E\ = 1 kcal/mol and a = 3.8 A. We design the L12 loop, disordered in the crystal 
structure, as a self-avoiding chain, employing the same Hamiltonian as EqJ3] with Ay = for 
i = 240 — 255 and all other j, or vice versa, which makes the L12 loop neutral to other parts of 
kinesin monomer (see also the SI text). The coiled-coil interaction responsible for dimerization 
between two neck-helices (residues 341—370) is modeled using 



%=\ j=i 



<Y>>-> 

ij 



12 



ij 



ij 



JVb 



i=X j=l 



ij 



A Qa ), 



(4) 



with El a 

A act 
ij 



1.2 kcal/mol. We defined A° a = 1 only for u i,j > 341" with R°. < R c , otherwise 
0, so that the two motor domains other than neck-helix repel each other. Because 
the two monomers are identical as a fold, we should impose the "same" topological bias to 
the both heads, which is achieved by setting Ay(T) = Ay (5) and r° (T) = rfj(B) for all i 
and j. This condition was used when we have previously studied the role of internal strain on 
the regulation mechanism of kinesin dimer 19] . To study the stepping dynamics in particular, 
however, we want the neck-linker of the tethered head in a disordered state, so that the tethered 
head can search and reach the next MT-binding site. This is realized by making the neck-linker 
(residues 327-338) of the tethered head always repulsive to the residues constituting the neck- 
linker binding site (see Fig. IB). To generate the kinesin configurations with A = 1 and A = 0, 
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we retain (A = 1) or discard (A = 0) the neck-linker zipper contacts (contacts inside green circle 
in Fig. IB) of the MT-bound head. 

The interaction between the kinesin and the MT is designed using 



jv k n mt 

Vk-mt = 

i=l k=l 




ef MT U— J -X^f-) \ + l B k B T^e^ 



^k 



E E { £ ™ (4) + hk B T^e^ [ (1 - A* fc ). (5) 



=1 k=l 



Nut is the number of all the residues, belonging to the tubulin heterodimers, that are within 
the interaction range from kinesin when the MT-bound head is poised at the central tubulin 
(o) (see Fig. 1A). The topological information of the binding interface between the MT-bound 
head and the central tubulin (o) is replicated to other surrounding tubulins (a-h) using A* k , 
so that the tethered head feels the identical potential on every MT-binding site. The native 
contacts between the kinesin and MT are defined between the i, k pairs within R^' MT = 9 
A. A* k = 1 for the native pairs, and A* k = otherwise. In addition to the non-bonded in- 
teraction (e^ T = ef T = 1 kcal/mol), the electrostatic interaction between the kinesin 
and the MT is considered because of the large amount of net negative charge (—35 e) in 
each tubulin unit. To preserve the native contact distance as the one in the crystal struc- 
ture even with electrostatic potentials, we adjust the parameter Xik m the Eqj5] by choosing 
Xik = 2 + (l B k B T/6e%- MT )[l/r° k + l/l D ]e~ r ^ /lD . At r = r° k , the energy for the native pair 
with the opposite charges is given as V(r° k ) = —e^' MT — l B k B T (7/6r° k + l/l D ) e~ r ^^ lD . The 
strength of electrostatics is controlled by the salt concentration c, that determines the Debye 
screening length l^ = (87tZbc) -1 / 2 ~ (3/\/c) A where l B = 7 A and c is in the unit of M 
(mol/l). For a pair formed at r° k m 8 A, the energy stabilization due to the electrostatics is 
l B k B T (7/6r° k + l/l D )e~ T y iD w 0.14 kcal/mol at c « 1 M, but the same value becomes ~ 1.4 
kcal/mol at c ~ 0.1 M, which is comparable to Eh = 1-0 kcal/mol. At physiological condition, 
c ~ (0.1 — 0.2) M for monovalent salt. A careful consideration is, however, required in choosing 



the c near the MT surface because of the counterion condensation [55J. Near the highly charged 
rod, the counterion concentration can be much higher than that of bulk. Therefore, we choose 
c = 1 M, which makes the electrostatics negligible, throughout the whole simulations (see SI text 
with SI Fig. 10 for the details). Eqj5] accomodates nonspecific interactions due to electrostatics 
even for the non-native contacts. 

Sampling the free energy surface. We combine Monte Carlo (MC) and molecular 
dynamics (MD) simulations to efficiently sample the kinesin configurations on the MT. We first 



generate a ensemble of initial configurations by both pivoting [48| a random position of neck- 
linkers (residues 327-339 for the MT-bound head, residues 328-337 for the tethered head) and 
translating the center of position of the tethered head. The acceptance of a trial move is decided 
by standard Metropolis criteria with the potentials defined in Eqs. [3j HI and [5j The subsequent 
hybrid MC-MD simulations are performed from a host of initial kinesin configurations. During 
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the 1000 simulation steps, we perform the MC simulations for the first 50 steps, and integrate 
the subsequent 950 steps using a velo city- Ver let algorithm. The kinesin configuration is collected 
every 1000 step. 
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Supporting Information 



Preparation of kinesin and MT structure : The simulations of kinesin were performed, 
referenced to an ADP-complexed crystal structure of rat kinesin dimer (PDB code : 3kin) in 
which the neck-linkers of both monomers are in an ordered state. For the completeness, us- 
ing a self-avoiding chain, we filled the gap of the missing residues 240-255, whose sequence is 
sktgaegavld, corresponding to the L12 loop in the crystal structure. Unlike the monomeric 
kinesin, KIF1A, whose L12 loop contains many lysine residues implying an important role in 



the motility by interacting with negatively charged E-hook of the tubulin |49j, the L12 loop of 
rat kinesin has only one lysine. Presumably, the electrostatic between the L12 loop and MT 
surface is not as important as the one in KIF1A. As an initial configuration, one of the kinesin 
monomers was placed on a tubulin binding site of the MT, and the other monomer is tethered 
to the MT-bound monomer via coiled-coil association but away from the direct MT interaction 
range. The topological information of binding interface between kinesin and tubulin was ac- 



quired using the Hirokawa and coworkers' KIF1A and tubulin complex 50J (see also Ref. 19] 
for the detailed procedure). 

Computation of two-dimensional potential of mean force (PMF) : The multiple 
histogram reweighting technique 5l|, |52| was adopted to compute the two-dimensional PMF 
between the tethered kinesin head and MT at temperature T. For example, the 2-D PMF can 
be obtained at arbitrary values of T if the conformational states are well sampled over a range 
of T values. The probability of finding the kinesin head at position (x, z) at temperature T is 
given by 

V p--B/r ElkMg£if) 
P(x, z)(T) = Efc=1 ^ — (6) 

2sE,x,z e Y%=1 n k e^ F k- E )/ T k 

where K is the number of histograms, hk{E, x, z) is the number of states between (E,E + 5E), 
(x,x + 5x), and (z,z + 6z) in the fc-th histogram, rik = J2e x z hk(E, x, z), T k is the temperature 
in the simulations used to generate the k— th histogram. The free energy, F^, that is calculated 
self-consistently, satisfies 

-Fr/T r _ \^ -E/T r J2k=lhk(E,X, z) 



e 



C -E/T r Z^k=l n kK^^,^) /«n 

^ £f=i^- E)/Tfe ' 



The self-consistent equation for Z k = e _Ffe//Tfe converge to the final values of {Zk} starting from 
Zk = 1 [k = 1, 2, . . . , K). Using the hybrid MC/MD simulations, we sampled the conformational 
states over the range of 295 K< T < 356 K. Once P(x,z)(T) is obtained, the 2-D PMF as a 
function of (x, z) is given by 

AF(x,y)(T) = -k B TlogP(x,z)(T). (8) 

AF(x,y), AF(y,z), and AF(Qi nt ,Q p ) are similarly obtained, and the 1-D PMFs are easily 
reduced from the 2-D PMFs. 
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Estimate of the translational diffusion constant of kinesin motor domain : When 
performing the Brownian dynamics simulation, we decided the diffusion constant of monomer 
using Stokes-Einstein relation, 

D o = (9) 
birr/a 

where i] is the viscosity of water (»s lcP = 10~ 3 N/m 2 ■ sec), and a is the hydrodynamic radius of 
residue (a ~ 0.19 + 0.14 nm), so that the monomer diffusion constant at T = 300 K is D m rs 578 
fim 2 ' j 'sec = 5.8 x 10 _6 cm 2 / sec. Similar calculation for the spherical object of radius R ~ 4 nm 
(approximately the size of kinesin head) leads to Dk ~ 5.5 x 10~ 7 cm 2 j sec = 55 fim 2 /nsec. 

Brownian dynamics with hydrodynamics : To simulate the real time kinetics of ki- 
nesin's swiveling motion using a Brownian dynamics of coarse-grained model, the inclusion of 
hydrodynamics is extremely important to obtain a correct time scale for the translational diffu- 
sion of whole object. Without hydrodynamics, the translational diffusion constant of a protein 
{Dk) that is coarse-grained by the N beads is scaled as D a /N, where D Q is the diffusion con- 
stant of the single bead. To naturally satisfy the Stokes-Einstein relation for the translational 
diffusion constant, the Brownian dynamics simulations requires the inclusion of hydrodynamic 
interaction. When off-diagonal elements of diffusional tensor is included and preaveraged, the 



diffusion constant of whole object Dk scales as 53|, [54 



Dk = ± + J^„± + £±^El. (io) 

N Qtx7]R h N N» N» v ' 

Dk ~ D /N u is the correct scaling relative to the D . Without hydrodynamics, the translational 
diffusion constant is significantly underestimated especially when N is large. However, it is also 
challenging to simulate the Brownian dynamics with hydrodynamic interaction. The Langevin 
equation under multidimensional space is 

n(t + At) - Ti (t) = j^fJ2 B v ■ F i + ^E B ^ • n iW ( n ) 

j j 

with (rij(i)) = 0, (ih(t) ■ rij(t)) = 5ij5(t — t') and D = BB T . The diffusion tensor D should be 
positive definite i.e. ■ Dy • Fj > for all F ^ 0, then D is deomposed into B and B T 

using the Cholesky decomposition. For the hydrodynamic diffusion tensor, Rotne-Prager diffu- 
sion tensor 54j is used. By averaging the square displacement of many different trajectories, we 
show in Figj9] how the explicit inclusion of hydrodynamics can alter the translational diffusion 
constant of the kinesin head domain. The value of D(hydro) = 19.9 fim 2 j 'sec is closer to the 
expected value for the -Dk(~ 55 /im 2 ' j 'sec. See above) that can be roughly estimated from the 
Stokes-Einstein relation. 

Consideration of electrostatics — Manning counterion condensation around the 
microtubule : The microtubule is viewed as a cylindrical object with many charges. Unlike 
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other thin polyelectrolyte such as RNA, ss-DNA and ds-DNA that can adapt its conformation 
depending on the salt concentration, the MT can serve as an excellent example to apply the 
Manning condensation theory for the cylindrical object. Because each tubulin hetero-dimer 
contains ~ 35e negative charge, the line charge density of microtubule is computed as 

d „(13x35e) A 

80A ' y ' 



The large Manning condensation parameter 55] £ = l B /b(~ 7.lA/(l/5.6)A) ~ 39.8 suggests 



that a drastic counterion condensation should occur around the microtubule to make £ ~ 1. 
This requires n ~ 34. le monovalent positive counterion, which is estimated from b* ~ Ib, i-e., 



: „ , 8 2 — r ~ 7.1, should condense to the surface of tubulin. The numerical solution of nonlinear 

13x (35— n) 

Poisson-Boltzmann equation with (j)'(a) = — — and 4>(oo) = 

V 2 = K 2 sinh0, (13) 



where 0(r) = ei/j(r) /ksT and k = vo7tTbc, determines the ion distribution around the charged 
cylinder. With 0(r) value, the number of positive and negative ions per around the cylinder 
can be computed using n + (r) = n^e'^ and n_(r) = n^e^ , respectively. The net condensate 
charge number q(r) = (n + (r) — n_(r)) provides the thickness of condensate ion by imposing the 
Manning condensation condition, q{Ru) x Ib = 1. For the parameter a = 12.5 nm, c = 150 
mM, the counterion condensation occurs at Rm ~ 16.7 nm (see Fig JTQj) . About 4 nm con- 
densate counterion layer is formed around the cylinder. Unless the interfacial binding takes 
place and expulges the counterion condensated near the binding site, the electrostatic potential 
due to the microtubule charge does not significantly affect the swiveling dynamics of kinesin 
head. The widely-accepted notion of "electrostatic steering" in the context of protein-protein 
association dynamics is not clear for the object on highly charged but effectively neutralized 
cylindrical surface. It has been shown that the processivity of the kinesin is affected by the salt 



concentration 
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Figure 1: Topology of MT, kinesin, and their interface. A. Geometry of MT surface showing the 
distances between the neighboring binding sites. \6a\ = \6e\ = 8.0 nm, \ob\ = \of \ = 10.9 nm, 
\6c\ = \dg\ = 6.3 nm, \od\ = \oh\ = 9.4 nm. The distances are measured using the position of the 
residue 400 at each a-tubulin subunit. B. The native contact map of kinesin model. The native 
bias between the neck-linker and its neck-linker binding site (neck-linker zipper contacts inside 
the green circles) is retained for the MT-bound head to have the neck-linker ordered (A = 1) while 
it is removed for the tethered head to make the neck-linker of the tethered head disordered. For 
structures with a disordered neck- linker in the MT-bound head (A = 0), the neck-linker zipper 
contacts are made repulsive. The contacts for the coiled-coil are colored purple. C. Native 
contacts between the residues of kinesin (yellow) and MT-binding sites (magenta). 
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Figure 2: PMF between the tethered head of kinesin and the MT surface for A = 1. B. 1-D 
PMFs in x-, y-, z- projections are shown at various temperatures. C. 2-D PMFs in xz-, yz-, xy- 
projections at T = 295 K. The binding sites e, d, and entropic basin S are marked with circles 
and arrows. D. 2-D PMFs at T = 355 K. 




Figure 3: PMF between the tethered head of kinesin and the MT surface for A = 0. B. 2-D PMFs 
in xz-, yz-, xy- projections at T = 295 K. The binding site c, which results in an intermediate 
state, is marked with circles. C. 1-D free energy profile F (x) is fitted by a harmonic potential 
1/2 kx 2 where k ~ 0.02 /c^T/nm 2 (dashed line). 
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Figure 4: The conformation of kinesin dimer (view from the back (A), top (B), and front (C)) 
when the tethered head is trapped at binding site c. The tethered head are partially bound to 
the site c, using only 1/3 of the native contacts. 




Figure 5: Results of the Brownian dynamics simulations of quasi-particle representing the teth- 
ered head on the PMF define in EqJTJ A. Conceptual representation of kinesin's stepping dy- 
namics (see text). B. The fraction of trajectories (f direct) directly reaching the target binding site 
e as a function of t p . C. (Left) The ensemble average of 100 trajectories generated for varying r p . 
Substeps are manifested for (x)(t)s with r p > 20 fis. The inset shows the substep in the (x)(t) for 
t p = 100 /is. (Right) An actual time trace Xi(t) is plotted in black with the time interval of 0.1 
jj,s from which the ensemble average is obtained as (x)(t) = X^=i x i(t)- The ensemble average of 
trajectories for r p = 20.0 /is is fitted to (x)(t) = 14.7 nmx [1 - 0.74e"' /15 - 3 ^ - 0.26e"* /149 ^] - 6.0 
nm. 
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Figure 6: PMF between the tethered head and the MT-binding site e as a function of two order 
parameters, Q p and Qi n t- A. (Top) The MT-binding motifs of the kinesin are colored in orange. 
(Bottom) An exemplary binding trajectory as functions of Q p and Qint- B. AF(Q int ), 1-D 
free energy profile as a function of Q int at varying temperatures. C. AF(Q int ,Q p ) at T = 305 
K(= T r ). As the binding progresses, minor structural disruption is observed. The free energy 
difference is color-coded from blue to red in k,BT r unit. D. AF(Q int , Q p ) at T = 355 K. The 
partial unfolding along the binding process is more pronounced than at a lower temperature. 
The overall binding pathway is drawn with a curved arrow. 
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Figure 7: Mechanochemical cycle of the conventional kinesin at the track of MT protofilament. 
ATP binding to the leading head (red) induces a undocked-to-docked transition of the neck- 
linker, which results in 16-nm stepping dynamics of the trailing head (blue). When the tethered 
head succeeds in binding to the next binding site, the rearward tension built on the neck-linker 
perturbs the nucleotide binding site and eases the dissociation of ADP from the tethered head 



571 . Il9| . The ATP hydrolysis follows at the trailing head (red). The step corresponding to the 



stepping motion is enclosed by a box. 
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Figure 8: A snapshot of simulation shown in the xz, xy, yz plane. The kinesin monomer in red 
is the MT-bound head, and the kinesin monomer in blue is the tethered head. 




t [rose] t [nsec] 

Figure 9: The diffusion constant of kinesin head domain is computed using Brownian dynamics 
simulations with explicit hydrodynamic tensor (Rotne-Prager tensor). Mean square displace- 
ments ((5R 2 )) of single bead (left) and kinesin head composed of 334 beads (right) are compared. 
For kinesin head, we performed Brownian dynamics simulations with and without hydrodynam- 
ics. By including the hydrodynamics, the correct diffusion behavior of kinesin head is obtained. 
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Figure 10: The numerical solution of nonlinear Poisson-Boltzmann equation under varying con- 
centrations of monovalent salt. The electrostatic potential around cylinderical geometry <j)(r) 
is obtained by solving V 2 = k 2 sinh with boundary condition, and 0(r) is used to compute 
n + (r), n_(r), and q(r). The inset and the purple line are the numerical solution for c = 150 
niM. The thickness of condense ion layer, Rm, is determined using the Manning condensation 
criterion, q(Rm) x l B = 1. 
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